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Abstract This paper discusses the dimensions of the spline spaces over T-meshes with lower degree. 
Two new concepts are proposed: extension of T-meshes and spline spaces with homogeneous boundary 
conditions. In the dimension analysis, the key strategy is linear space embedding with the operator of 
mixed partial derivative. The dimension of the original space equals the difference between the dimension 
of the image space and the rank of the constraints which ensuring any element in the image space has a 
preimage in the original space. Then the dimension formula and basis function construction of bilinear 
spline spaces of smoothness order zero over T-meshes are discussed in detail, and a dimension lower bound 
of biquadratic spline spaces over general T-meshes is provided. Furthermore, using level structure of 
hierarchical T-meshes, a dimension formula of biquadratic spline space over hierarchical T-meshes are 
proved. A topological explantation of the dimension formula is shown as well. 

1 Introduction 

A T-mesh is a rectangular grid that allows T-junctions. T-splines, a type of point-based splines defined over 
a T-mesh, were proposed by T. W. Sederberg in JHJ [9J, and have become an important tool in geometric 
modeling and surface reconstruction. For T-splines, T-mesh plays two roles: defining its parametric domain 
decomposition fashion and representing the topological structure of the control net of a T-spline surface. 

According to its definition, a T-spline is a piecewise polynomial, instead of a single polynomial, within 
a cell of the T-mesh. This is incompatible with the standard defining fashion of classic splines. Recall that, 
given a spline knot sequence, a univariate spline can be defined, which is a single polynomial between any 
two neighboring knots. Hence in some of the present authors with some other coauthors introduced the 
concept of spline spaces over T-meshes, where every function in the space is exactly a polynomial within 
each cell of the T-mesh. A dimension formula can be proved with the B-net method [3 1 and the smoothing 
cofactor method [5 1 for the spline space S(m, n, a, (3, 3?) as m ^ 2a + 1 and n ^ 2/3 + 1. Then in H we 
provided an approach to define the basis functions of the bicubic splines with smoothness order one. The 
applications of the basis functions in surface fitting is explored as well. 

According the dimension formula in J3), one can obtain the specified dimension formulae for some 
spline spaces with low degree as S(l, 1, 0, 0, ST), S(2, 2, 0, 0, ST), S(3, 3, 0, 0, and S(3, 3, 1, 1, 
Furthermore, with similar approaches in (HE), one can construct their basis functions with some "good" 
properties, say, compact support, nonnegativity, and forming a partition of unity. In order to achieve high 
order smoothness with as low as possible degree, we expect to obtain the dimension formulae and basis 
functions construction for the spline spaces S(m,n,m — l,n — 1, Especially, the most interesting 
ones are those for S(2, 2, 1, 1, 5?) and S(3, 3, 2, 2, In the paper, we only focus on the dimension of the 
former space. 

In the following, we first introduce two concepts: the spline space with homogenous boundary con- 
ditions (HBC) and the extended T-meshes associated with some spline space. As a foundation of the 
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later analysis, we discuss in detail the dimension formula and basis functions construction for the space 
S(l, 1, 0, 0, 2?Y In J3) we have shown that the dimension of a bilinear spline space is the sum of the num- 
bers of crossing vertices and boundary vertices in the given T-mesh. However, the proof proposed here 
shares the same fashion as in the analysis of a lower bound of dimensions of biquadratic spline spaces, and 
avoids the problem of recycling dependence of T-vertices. 

An important technique in the dimension analysis is linear space embedding by an operator of mixed 
partial derivative, which embeds the space S(m, n, m— l,n— 1, 2F) into the space S(m— 1, n— 1, m— 2, n— 
2, A necessary and sufficient condition for describing any element in S(m— 1, n — 1, m — 2, n — 2, 
is the image of an element in S(m, n, m — l,n — 1, by the operator. In the paper, we just discuss cases 
of m = n = 1 or 2. With the method, a lower bound of dimensions of S(2, 2, 1, 1, 3F) is proved. Finally, by 
making use of the level structure of hierarchical T-meshes, a dimension formula is provided for biquadratic 
spline spaces over hierarchical T-meshes. 

The paper is organized as follows. In Section 2, the spline spaces over T-meshes is reviewed and two 
concepts are proposed: extension of T-meshes and spline spaces with homogeneous boundary conditions. 
Then the dimension formula and basis function construction of bilinear spline spaces of smoothness order 
zero over T-meshes are discussed in detail in Section 3, and a dimension lower bound of biquadratic spline 
spaces over general T-meshes is provided in Section. In Section 5, using level structure of hierarchical 
T-meshes, a dimension formula of biquadratic spline space over hierarchical T-meshes are proved. A topo- 
logical explantation of the dimension formula is shown as well. Section 6 concludes the paper with some 
discussions. 



2 T-meshes and Spline Spaces 

A T-mesh is basically a rectangular grid that allows T-junctions. As for details of T-meshes, please refer to 
0. Here the T-meshes we discuss are regular, i.e., the domaines occupied by the T-meshes are rectangles. 
We adopt the compatible definitions of vertex, edge, cell with those in |3|. 
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Figure 1 : Horizontal T-vertices (a and b) and vertical T-vertices (c and d) 



T-vertices in a T-mesh can be classified into two types: horizontal T-vertices and vertical T-vertices. 

The T-vertices as shown in Figures Q]a and b are horizontal ones, and those as shown in Figures[T]c and d 
are vertical ones. In [3 1, edges and c-edges are defined. Here we introduce a new type of edges, 1-edges. A 
horizontal/vertical 1-edge is a continuous line segment which consists of some horizontal/vertical interior 
edges and whose two endpoints are boundary vertices or horizontal/vertical T-vertices. In other words, an 
1-edge is a longest possible line segment in the mesh. The boundary of a regular T-mesh consists of four 
1-edges, which are called boundary 1-edges. The other 1-edges in the T-mesh are called interior 1-edges. For 
example, in Figure [2] the given T-mesh 3? has three horizontal interior 1-edges and three vertical interior 
1-edges, respectively. 

Given two series of real numbers Xi, i = 1, . . . , m, and yj, j — 1, . . . , n, where x% < Xi+i and 
Uj < yj+i, a rectangular grid can be formed with vertices (xi,yj), i — 1, . . . , m, j — 1, . . . , n. This 
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grid is called a tensor-product mesh, denoted by (x\, . . . , x m ) x (yi, . . . , y n ), which is a special type of 
T-mesh. From a regular T-mesh ST, a tensor-product mesh 5 TC can be constructed by extending all the 
interior 1-edges to the boundary. 2? c is called the associated tensor-product mesh with 3*. See Figure [2] 
for an example. 
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Figure 2: A T-mesh with its associated tensor-product mesh S7 C and its extension 3 re . 



2.1 Spline spaces over T- meshes 

Given a T-mesh 3T, we use & to denote all the cells in 2T and fl to denote the region occupied by all the 
cells in & . In [3|, the following spline space definition is proposed: 

S(m,n,a,0,&) := {f(x,y) g C a > p (n) : f(x,y)\+ g P„ m ,V0 g (2.1) 

where P m „ is the space of all the polynomials with bi-degree (m, n), and C a,/3 is the space consisting of 
all the bivariate functions which are continuous in £1 with order a along x direction and with order [3 along 
y direction. It is obvious that S(m, n, a, (3, 3?) is a linear space. 

Now we introduce a new spline space over T-meshes with the following defintion: 
S(m, n, a, (3, ST) := {f(x, y) g C a,/3 (R 2 ) : f(x, y)\* g P mn , V0 g J?, and /| Ra \ n = 0}, (2.2) 
which is called a spline space over the given T-mesh & with homogeneous boundary conditions (HBC). 

2.2 Extended T-meshes 

In order to discuss the dimension of the spline space S(m, n, m — 1, n — 1, 3?), we extend the given regular 
T-mesh 3~ in the following fashion: Pickup a tensor-product mesh jtfi with 2(m + 1) vertical lines and 
2(n + 1) horizontal lines, such that the central rectangle of j$ is the same size as SI, the region occupied by 
3*"; Then all the edges with an endpoint on the boundary of & are extended to reach the boundary of M '. 
The result mesh, denoted as £? e , is called an extension of the original T-mesh (or say, an extended T-mesh) 
associated with the present spline space. Figure |2] provides an example, where the T-mesh is extended 
associated with S(2, 2, 1, 1, 3 7 ). 

It should be noted that the extension is associated with a spline space. Hence, associated with different 
spline spaces, one will obtain different extended T-meshes. The following theorem shows that, using an 
extension of T-mesh, the dimension analysis of S(2, 2, 1, 1, 3? e ) is the same as the dimension analysis of 
S(2,2,l,l,^). 
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Theorem 2.1 Given a T-mesh 3?, assume S? z is its extension associated with S(2, 2,1,1, Then 

S(2, 2, 1, 1, &) = S(2, 2, 1, 1, (2.3) 
dimS(2,2,l,l,^) = dim S(2, 2, 1, 1, ST e ). (2.4) 

Proof: We first prove Equation (12.31 >. Since 3~ e is an extension of ST, it follows that 

S(2,2,l,l,^ e )|^ c S(2, 2,1,1,^). 

In the following we will prove that, for any / € S(2, 2, 1, 1, 5"), there exists / g S(2, 2, 1, 1, £f e ) such 
that/|^=/. 
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Figure 3: Bezier ordinates in an extended T-mesh 

Here we define the function / over ST e by assigning its Bezier ordinates in every cell of S? £ \,5?. Since 
/ e C 1,1 , and there does not exist any interior T- vertices in ,S? £ \£?, f shares the same ordinates along 
the common boundary between two neighboring cells. (Therefore, these same ordinates appears once in 
Figure [3] which illustrates the distribution of the Bezier ordinates in the left-bottom part of the extended 
region £/' e \£/ r .) According to the specification of / over 3 ', it follows that the ordinates labeled with 
"•" are determined. Since / meets the zero function along the boundary of the extended T-mesh S7 e with 
order one, the ordinates labeled with "A" are determined as well. The rest ordinates labeled with "□" can 
be determined by its corresponding neighboring horizontal or vertical ordinates labeled with "A" and 
Here the continuous condition is the collinear of the corresponding three points. It should be noted that the 
ordinates themselves determined in this fashion ensures that they define a continuous spline automatically. 

After determining the former ordinates, the ordinate labeled with "■" can be determined as follows: It 
follows that the eight ordinates neighboring the ordinate "■" define a bilinear function. When the ordinate 
corresponding to "■" also lies on the bilinear function, the continuous of order one along two directions 
can be guaranteed. In fact, when dealing with the determination of all the ordinates in the extension region 
around ST, if we meet cases of recycling determinations, we can address the ordinates using this method of 
bilinear functions construction. 

Hence we obtain the specification of the function / in S(2, 2, 1, 1, 3? e ), which satisfies f\& = f. 
According to the former analysis, it follows that 

S(2,2,1,1,^) = S(2,2,1,1,^ £ )|^. 

In order to prove Equation (12.4b . we will prove that, for any nonzero function / 6 S(2, 2, 1, 1, £? e ), it 
follows that f\g- ^ 0. Suppose =0, i.e., there exists p = (xo, yo) 6 ^ 6 \^ such that f{xo,yo) 0- 
And assume ft = (x u x r ) x (y b ,y t )- If xi ^ x ^ x r , y < y b , then f(x ,y), y ^ y b is a quadratic 
spline, whose support has just three breakpoints. This contradicts with the fact about the minimal support 
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of quadratic splines (The support of a nonzero quadratic spline should have at least four breakpoints). 
Similarly, we can show that p is outside of the regions xi ^ xq ^ x r , yo > yt\ yi ^ yo ^ y r , %o < %i an d 
yb ^ yo ^ Vu a^o > x r- Therefore, without loss of generality, we assume xq < x r , yo < yb- But this also 
results in a nonzero quadratic spline with three breakpoints in its support. Hence we obtain f\& ^ 0, and 
has proved that 

dim S(m, n, a, j3, 3?) = dim S(m, n, a, /?, 3 re ). 
This completes the proof of the theorem. I 

According to Theorem 12.11 we will consider only the biquadratic spline spaces over T-meshes with 
HBC. 

Remark: A similar analysis can show that, for any m and n, it follows that 

S(m, n, m — 1, n — 1, 2?) = S(m, n,m — l,n — 1, 5" e )|5>, (2.5) 
dim S(m, n,m — 1, n — 1, 5") = dim S(m, n, to — 1, n — 1, 5" £ ), (2.6) 

where ^ e is an extension of 5^ associated with S(m, n, to — 1, n — 1, ^7). Hence one can replace the 
dimension discussion for S(m, n,m — 1, n — 1, ^) with that for S(m, n,m — l,n — 1, 



3 Dimensions and Basis Functions of S(l, 1, 0, 0, ^) 

In J5] we have proved the following dimension formula of the spline space S(m, n, a, f3, 3") over T-mesh: 

dimS(m, n, a, (3, ST) = F(m + l)(n + 1) - E h (m + l)(/3 + 1) 

-E v [n + l)(a + 1) + V(a + 1)09+ 1), (3.1) 

where m 5^ 2a + 1, n ^ 2/3 + 1, F is the number of all the cells in 2T, Eh and E v the numbers of horizontal 
and vertical interior edges, respectively, and V the number of interior vertices (including crossing and T- 
vertices). Specially, as m = n = 1 and a = (3 = 0, it follows that 

dimS(l, 1,0,0, ST) = V+ + V b , (3.2) 

where V + is the number of crossing vertices, and V b the number of boundary vertices. Using an extension 
of the T-meshes associated with S(l, 1,0,0, 3*), it follows that 

dimS(l,l,0,0,^ e ) = dim S(l, 1,0,0,5") = V+ , (3.3) 

where is the number of the crossing vertices in the extended T-mesh £? £ , which equals the sum of the 
numbers of crossing vertices and boundary vertices in 3^. 

Now we will prove that the former formula holds for a general regular T-mesh (not just an extended 
T-mesh). The method taken in the proof is similar with the method proposed in the next section for proving 
the dimension properties of S(2, 2, 1, 1, 3*). This method solves the problem of recycling dependence of 
T- vertices which happens in the proof proposed in [3 1. 

Now we introduce some notation and lemmas for the proof. In the given T-mesh 3?, let E denote the 
number of interior 1-edges, and V + the number of crossing vertices. Define S(0, 0, — 1, — 1, 3?) to be the 
space consisting of functions which are a constant in each cell of the T-mesh 3?, and have no smoothness 
requirement between neighboring cells. It is obvious that its dimension is the number of all the cells in the 
T-mesh. 
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3.1 Operator of mixed partial derivative 



The operator of mixed partial derivative is introduced as follows: 
d 2 - 

V := — — : S(m. n,m-l,n- 1, J) -> S(m - 1, n - 1, m - 2, n - 2, (3.4) 
oxoy 

where m, n ^ 1. Since the function satisfies HBC, the operator is a one-to-one, but not onto mapping. For 
example, for any nonzero g G S(to, n, m — 1, n — 1, ^"), it follows that £>(<?) must reach positive values in 
some parts and negative values in some other parts. Hence a nonnegative function in S(m — 1, n — 1, m — 
2, n — 2, ^) has no pre-image under the operator X>. 

Define ^ 

l(g)(x,y) = [ [ g(s,t)dsdt. (3.5) 



-co J — oo 



It follows that, for any / G S(m, n, m — 1, n — 1, 5^), I o P(/) = /. Hence T can be considered to be the 
inverse operator of £>. 

It follows that, for any g G S(m — 1, n — 1, m — 2, n — 2, 1(g) is a piecewise polynomial of degree 
(m, n) with smoothness C""- 1 -."- 1 , With respect to functions in S(m, n,m — l,n — l, d m T(g)/dx m 
or d n I(g) / dy n may be discontinuous inside some cells of T. But the discontinuity must happen on the 
extension of some 1-edges of T. Let 

d m ,n = dim S(m, n,m — l,n — 1, i^). 

If we have known d m _i „_i, i.e., the dimension of S(m — l,n — 1, m — 2, n — 2, ^T), and the number 
r m -i, n -i of linear-independent constraints ensuring G S(m,n,m — l,n— 1, ^) for any g G S(to — 
1, n — 1, m — 2, n — 2, ^"), then it follows that 

dm,n — dfxi— 1,ti — 1 Tin— l.n— 1- (3.6) 

On the other hand, if there are r' m _ x n _ x constraints, which may be linear dependent, proposed for ensuring 

1(g) G S(m, n, to — 1, n — 1, ^) for any g G S(m — 1, n — 1, m — 2, n — 2, c^ 7 ), then it follows that 

dm,n ^ d m — x,n— 1 ^m— l,n— 1' (3-7) 

Equations (13.61 > and (13.7b are used to prove the dimension formulae of bilinear and biquadratic spline spaces 
over T-meshes in the rest of the paper. 

3.2 Some lemmas 

As m = n = 1, the following lemma proposes the constraints ensuring 1(g) G S(l, 1, 0, 0, 3F) for any 
g eS(0, 0,-1,-1,^). 

Lemma 3.1 Given a regular T-mesh ST, let g G S(0, 0, —1,-1, J 7 ). 77zen 

T( 9 )(x,y)eS(l,l,0,0,J) 
if and only if the following two sets of conditions are satisfied simultaneously: 
1. For any horizontal l-edge l h , 

g(s,y -)ds = g(s,y +)ds, (3.8) 
where the two end-points ofl h are with coordinates (xq, yo) and (x±, yo); 
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2. For any vertical l-edge V, 



Vi rVi 

g(x -,t)dt= g(x +,t)dt, 

Vo J Vo 



(3.9) 



where the two end-points ofl v are with coordinates (xq, yo) and (xq, yi). 



Proof: For any horizontal l-edge l h , its right end-point (x\, yo) is a T-vertex or a boundary vertex. Suppose 
the edge through the vertex (xx, Vo) and perpendicular with l h is e v (as shown in Figure @). If 2(g) (x, y) £ 
S(l, 1, 0, 0, ,'7), then 2(g)\ e v is a linear polynomial in a neighborhood of (xx, yo)- Hence it follows that 



l" 
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- (xi,yo) 



Figure 4: A horizontal l-edge (1) 



d d 
—X(g){x 1 ,y -) = —l(g)(xi,y +). 

According to the definition of 2 in Equation (13.5b . Equation ( 13.101 can be rewritten into 



j(s,y -)ds 



g(s,y Q +)ds. 



(3.10) 



(3.11) 



Extend l h to reach the left boundary of the T-mesh. If there are not other horizontal 1-edges on the extension 
edge, then, in every cell that intersects the extension edge, g is a pure linear polynomial. Then g(s, yo — ) = 
g(s,y +) as s < x , and 

g(s,y -)ds = / g(s,y +)ds. (3.12) 



According to Equations ( 13.1 U and (13.121 i. it follows that the equation (13.8b in the first condition holds for 
any horizontal l-edge. 

It there exist some other horizontal 1-edges on the extension edge, then we run through these horizontal 
1-edges from left to right. For the first one, we can prove that Equation (13 - 8b holds with the former approach. 
Then we consider the other 1-edges one by one, and prove that the corresponding equation ( 13.121 ) holds as 
well, and then Equation (13.8b holds. Hence we finish to prove Equation (13.8b holds for every horizontal 
l-edge. 

For any vertical l-edge, we can prove that the corresponding equation ( 13. 9b holds in a similar way. 
Hence the necessity of the lemma is proved. 

Now we will prove the sufficiency of the lemma. Suppose g G S(0, 0, — 1, — 1, 3F) satisfies the two 
sets of conditions in the lemma, but 1(g) ^ S(l, 1,0,0, 3T), i.e., there exists a cell where dl(g) / dx or 
dl(g)/dy has at least a discontinuous point. Define two sets consisting of the cells of & as follows. Let 
B x denote all the cells where 31(g) / dx has at least a discontinuous point; Let B y denote all the cells where 
82(g) jdy has at least a discontinuous point. Then according to the assumption, B x U B y is not empty. 
Without loss of generality, we assume B y is not empty. Consider a cell in B y whose left-bottom corner 
has the minimal y coordinate in B y . If there exist more than one such cell, select one from them such that 
its left-bottom corner has the minimal x coordinate. Then now we have selected a unique cell c. In c, 
dl(g) I dy has at least one discontinuous point (x, y). Hence it follows that 



g(s,y-)ds ± \ g(s 7 y+)ds 



(3.13) 
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Consider the horizontal straight line through (x,y). On the straight line if there does not exist an 1-edge of 
2? on the left side of (x, y), then Equation ( 13.13b fails to hold. If there exists at least one 1-edge on the left 
side of (x, y), then according to the first condition, a contradiction arises as well. Hence B y is empty. In the 
same fashion, it follows that B x is empty as well. This contradicts with B x U B y non-empty. Therefore we 
finish the proof of the sufficiency. I 

There are E + 4 conditions in Lemma [3~T1 The following lemma states that these E + 4 conditions are 
equivalent with E + 2 conditions in another form. 

Lemma 3.2 Given a regular T-mesh the occupied rectangle by ST is (xi,x r ) x {yb,yt)> where the 
different y coordinates of horizontal l-edges are y < yy < . . . < y n . Suppose g £ S(0, 0, — 1, — 1, 5 r \ 
Then 

/>X r />X r 

/ 9{s,Vi-)ds= g(s,yi+)ds, i = Q,l,...,n (3.14) 

Jx l Jxi 

is equivalent with 

g(s,y)ds = 0, y e j/i+i), i = 0, . . . , n - 1. (3.15) 



It should be noted that, since g is a piecewise constant, J^ r g(s, y)ds, y S (yi, yi+i), is a constant inde- 
pendent on y. 

A similar statement can be made for vertical l-edges. 
Proof: Since g is a piecewise constant in , and zero out of it follows that 

g(s, y -)ds = / g(s, y)ds = 0, y < y , 

J XI 

rx r px r 

/ g(s,yo+)ds= g(s,y)ds, y <y<Vi- 

Jxi J Xi 

Hence according to Equation d3.14t as i = 0, we have 

g(s,y)ds = 0, y € (ya,Vi)- 



This means that Equation (13.15b holds as i = 0. Then we consider cases as i = 1, 2, . . . , n — 1. Similar 
discussions guarantee Equation ( 13.15b holds as i = 1, 2, . . . , n — 1. 

On the other hand, for any i = 0, 1, 2, . . . , n — 1, if all the quations in d3.15t hold, it is easy to prove 
that all the equations in (13.14b hold. H 

Recall that, for a given T-mesh 3^, an associated tensor product mesh 3^ c can be obtained by extending 
all the interior l-edges to the boundary of 3T. Suppose there are E 1 interior l-edges in 3~ c and E interior 
l-edges in ST , where E' ^ E (since it is possible that more than one l-edges in ST on the same 1-edge in 
3 rc '). Select any horizontal 1-edge I from 3~ c , and suppose the horizontal l-edges ?i, . . . , Ik in 3^ lie on I. 
Then it follows that the constraints 



g(s,y-)ds = / g(s,y+)ds, i = 1, 
is equivalent with the constraints 

g(s,y-)ds= / g{s,y+)ds, i = 1, . . . , fc — 1, I g(s,y-)ds = I g(s,y+)ds, 



where y is the vertical coordinate of I. According to Lemma l3T2l all the following constraints, with respect 
to any horizontal 1-edge I in 3^ c , 

g(s,y-)ds = / g(s,y+)ds 
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are equivalent with the integrations of g along all the span formed by two neighboring horizontal 1-edges 
in ^ c are zero. A similar equivalence can be made for vertical 1-edges. Here we should noted that, for 
any boundary 1-edge, there is only one constraint stated in Lemma [3T| Hence the necessary and sufficient 
conditions that 1(g) E S(l, 1, 0, 0, &) for g 6 S(0, 0, -1, -1, ST) are E + 2 constraints in this fashion. 
On the other hand, if = (xq, x\, . . . , x m ) x (yo, Vi> ■ ■ ■ , Un)> it follows that 



I g(s, t)dsdt = (x i+ i - x l )C l = ^T(yj+i - Vj)Dj, 

-oo J —oo : n A n 



m — 1 n — 1 



where 



i=0 j=0 



Ci= [ g(x,y)dy, x £ (a*, x i+1 ), 
J yo 



Dj = / g{x,y)dx, y E (y 3 ,y j+1 ). 
Jxo 

Hence these E + 2 constraints is with defective rank at least one. The latter dimension theorem [3~4l will 
show that the defective rank is exactly one. 

Up to now, we have finished the description of the constraints that 1(g) E S(l, 1, 0, 0, ST). Before we 
state the dimension theorem, a topological equation of T-meshes is proposed in the following lemma. 



Lemma 3.3 Given a regular T-mesh 3T, suppose 3? has F cells, V + crossing vertices, and E interior 
l-edges. Then 

F = V + + E + 1. 

Proof: Suppose, in ST, there are V T T-vertices, V bT boundary vertices (excluding four corner points). 
Since every cell has four vertices, running through all the cells will meet every crossing vertices four times, 
every interior T-vertices twice, every corner points once, and every other boundary vertices twice. Hence it 
follows that 

4F = W + + 2V T + 2V bT + 4. 

On the other hand, the end-points of every interior 1-edge are either interior T-vertices or boundary vertices 
(not corner points). Therefore, we have V T + V bT — 2E. From these two equations one gets F = 
V+ + E + 1. ■ 

3.3 Dimension theorem 

Theorem 3.4 Given a regular T-mesh ST with V + crossing vertices, it follows that 

dim S(l, 1,0,0,^) = V + . 

Proof: We first prove that dim S(l, 1, 0, 0, 3~) V~^~ . Suppose a function / E S(l, 1, 0, 0, 3?) reaches 
zero at all the crossing vertices. We intend to prove that / = 0. 

Suppose / ^ 0. Without loss of generality, we assume / is greater than zero in some regions in ST . 
Then there exists a point p in such that f(p) = S = max / > 0. Let p be in the cell c. There will only 
happen the following two cases: 

1 . there exists an edge e of c such that / is a constant S along e; 

2. for any edge e of c, / is not a constant along e. 
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Consider Case 1. Let 1$ denote the 1-edge on which e lies. Define a set L which consists of all the 
1-edges on which / are constants 5. P consists of all the end-points of the 1-edges in L. Since Iq 6 L, 
both L and P are non-empty. Now we category the vertices in P into two types. If a vertex in P is also an 
interior vertex on another 1-edge in L, then the vertex is called a flat vertex. Otherwise it is called a non-flat 
vertex. In the following we will prove that there must exist at least one non-flat vertex in P. If so, select 
one non-flat vertex q. Then q must lie on an 1-edge which is not in L. q is an interior point on l\. Since 
f(q) = S and is not a constant, there exists a point r on l\ such that f(r) > S, which contracts with the 
fact that S is the maximum of / over S? . 

In fact, the non-flat vertex q can be selected to the vertex in P with the minimal y coordinate. If there 
exist more than one vertex with the minimal y coordinate, one with the minimal x coordinate among them 
is selected. Such the selection ensures that q is a non-flat vertex. If not, q lies on two 1-edges 1% and Z3 in 
L, where l 2 is horizontal and l 3 vertical. If q is a horizontal T-vertex, then the bottom end-point of 1% has a 
smaller y coordinate than q. If q is a vertical T-vertex, then the left end-point of l 2 has the same y coordinate 
as q, but with smaller x coordinate than q. Both the cases contradict with the selection of q. Hence q is a 
non-flat point. Then we can prove / = for Case 1 . 

Now we consider Case 2. Since f\ c is a bilinear function and / is not a constant along any edges of 
c, / reaches its maximum only at one of its corners. Hence p is a corner of c, which is a T-vertex, say a 
horizontal T-vertex. The vertical 1-edge through p is assumed to be I4, which takes p as its interior point. 
Since f(p) = S and /|; 4 is not a constant, there exists a point s on I4 such that f(s) > S, which contracts 
with the fact that 5 is the maximum of / over 3? '. 

Summarizing the consideration for both the cases, it follows that / = 0. Then all the crossing vertices 
in & form a determining set of the spline space S(l, 1, 0, 0, According to the theory of the determining 
sets in spline functions JUG!, one gets 

dimS(l, 1,0,0, ST) < V + . (3.16) 



On the other hand, we will prove that dim S(l, 1, 0, 0, 8F) ^ V + in the following. To do so, we consider 
the operator of mixed partial derivative defined in Equation ( 13.41 i as m = n = 1: 

V: S(l, 1,0,0,^) ->S(0, 0,-1,-1,^). 

Here V is injective. The spline space S(0, 0, — 1, — 1, 3F) is with dimension F, the number of the cells in 
3T. According to the analysis in the former section, in order to ensure 1(g) G S(l, 1, 0, 0, one needs 
to satisfy E + 2 constraints, which have defective rank at least one. Then according to Lemma [331 

dim S(l, 1,0,0,^) ^ F - (E + 2) + l = V + . (3.17) 

Combining Equations ( 13.161 ) and ( 13.171 ). it follows that the dimension theorem is proved, i.e., 

dim S(l, 1,0,0,^) = V + . 



3.4 Basis functions 

We can construct a set of basis functions {bi(x, y)}Y=i f° r me spline space S(l, 1, 0, 0, ST). The basis 
functions should have the following properties: 

1. Compact Support: For any i, bi(x, y) has a support as small as possible; 

2. Nonnegativity: For any i, bi(x, y) ^ 0; 
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3. Partition of Unity: If 3" is an extension of some T-mesh ,% with respect to the spline space 
S(l, 1, 0, 0, 3), and the region occupied by ,% is SI, then 

y+ 

^2h(x,y) = 1, (x,y)eSl. 

i=l 

The basis functions with those properties can be constructed as follows: Suppose the crossing vertices in 
3 are Vi with coordinate (xi,yi), i = 1, . . . , y + . Then we require the function bi(x, y) satisfy bi(xj,yj) = 
5ij. According to the dimension theorem [3~4l and the first part in its proof, bi(x, y) is determined uniquely. 
All the functions bi(x, y) form a set of basis functions of the spline space, which have the former three 
properties. In fact, it is easy to show that Properties 1 and 2 are satisfied. Now we prove Property 3. 

Theorem 3.5 Given a regular T-mesh 3 ', which occupies a rectangle SI, its extension associated with 
S(l, 1, 0, 0, 3) is 3? e . The crossing vertices in S? e are Vi with coordinate (xj, j/j), i — 1, . . . , V + . The 
function set {bi(x, y)}Y=i C S(l, 1, 0, 0, 3 e ) satisfy bi(xj, yj) = Jy. Then 

v+ 

J2h(x,y) = l, (x,y)efl. (3.18) 

i=l 

Proof: Let 

v+ 

f(x,y) = ~J2h{x,y), 

i=l 

and I denote the boundary of 3. In order to show Equation (13.18b holds, we first prove f\i = 1. Because 
the vertices on £ are crossing vertices in 3 e , it follows that / reaches 1 on these vertices. Since f\e is a 
piecewise linear function with knots being these vertices, it follows that f\i = 1. Then, in the following, 
we prove that, for any (x, y) G CI, f(x, y) = 1. Consider the function 

/ x !f{x,y)-l fay) e SI 
g{x,y) — < 

I otherwise 

Since f\i = 1, one has g £ S(l, 1, 0, 0, 3"). But g is zero at all the crossing vertices of 3", it follows that, 
according to the proof of the dimension theorem [3~4l we have g = 0, i.e., f(x, y) = 1, (x, y) £ SI. I 

Remarks: There are some interesting problems open here. 

1 . How to directly specify the former basis functions in every cell of a general T-mesh? 

2. How to evaluate the function or the surface which is represented in the linear combination of the 
former basis functions? 

3. What is the "knot" insertion algorithm in this spline space? 

Though the space is just bilinear, the solutions to these problems will possibly hint us how to do in higher 
degree spline spaces over T-meshes. 



4 A Lower Bound of the Dimension of S (2, 2, 1, 1, 

We can apply a similar method proposed in the proof of the dimension theorem of S(l, 1, 0, 0, 3) in the 
former section to the dimension analysis of S(2, 2, 1, 1, 3). After that, we can obtain a lower bound of the 
dimension, i.e., 

dim S(2, 2,1,1,5') > V + -E+l. (4.1) 
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4.1 Some lemmas 



Now we consider the operator of mixed partial derivative as follows: 

d 2 



V 



dxdy 



S(2,2,l,l,^)-fS(l,l,0,0,^). 



Here T> is injective as well. The operator 1(g) is defined in the same way as Equation ( 13.5b . The following 
lemmas discuss the constraints ensuring 1(g) E S(2, 2, 1, 1, for any g £ S(l, 1, 0, 0, 3^). 

Lemma 4.1 Given a regular T-mesh 27 , let the coordinate of the end-points of any horizontal l-edges l\ be 
(x\ x , ]fl) and (x\ 2 , Ui), i = 0,1, ... ,m, and the coordinate of the end-points of any vertical l-edge ZJ be 
(xj, and (arj, yj 2 )> j = 0, 1, . . . , n. For any g 6 S(l, 1, 0, 0, 3T), it follows that 



I( s )eS(2,2,l,l,J)« 



<2 Q 

-rr9(s,y^-)ds 
rvh d 



,2 Q 

Q^g(s,y'i+)ds, £ = o, l, . 



, m, and 



fVj 3 Q fV 3 2 Q 

/. dx' g W-> t)dt= I v te9W+,t)dt, j = 0,l 
< y j i Vj i 



Proof: We first prove the necessity 



Let f(x,y) — I(g)(x,y). Without loss of generality, we only 



PMi.yl) 



Figure 5: A horizontal l-edge (2) 



prove that, when / 6 S(2, 2, 1, 1, 3^), the constraints corresponding to horizontal l-edges are satisfied. As 
shown in Figure|5] the two end-points of the horizontal l-edge are P\ (x^ , y\ ) and Pi (x\ 2 ,y\). The vertical 
edge on which P2 lies is /. Because / € S(2, 2, 1, 1, ST), it follows that f(x[ l 2 ,y) is a quadratic polynomial 
with respect to the variable y in a neighborhood of P 2 along I. Hence, 

/(4>^-) = /(4,^+). 



Q y f( x i2^yi- 
o 2 



dy- 



According to the definition of / and the continuous of g, the first two equations hold trivially. Substituting 
the definition of / into the last equation, one has 



« d_ 

x> dy 



9(s,Vi-)ds 



2 A 

-00 9y' 



,(s,y l l+)ds. 



(4.2) 



Extend the current l-edge to the left boundary of the T-mesh. If there does not exist any other l-edges 
on the extension, then in every cell which intersects the extension on the left side of (x^ x , y^), g is a single 
bilinear function. Hence it follows that 



d f Xil d 

—g(s,y?-)ds= / — g (s, $ +)ds. 
dy J-oo dy 



(4.3) 
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Subtracting Equation ( 14.31 ) from Equation ( 14.2b , one gets 

r x i2 g r x % g 

/ T r g{s,y?-)ds= / — g(s,y?+)ds. 

This proves that the equation in the lemma corresponding to the current 1-edge holds. If there exists some 
other 1-edges on the extension, then we consider these 1-edges one by one from left to right. We can prove 
that the former equations hold for these 1-edges. Hence for the current 1-edge, according to Equation (14.2b . 
the same equation holds as well. 

Now we prove the sufficiency <=. We will show that, for any g G S(l, 1, 0, 0, if the following 
equations hold: 

x * 2 d f x * 2 d 

—g(s,y?-)ds = / —g(s,y?+)ds, i = 0,1, •••,"», (4.4) 
h dy dy 

vh g rvh g 

—g(x--,t)dt = / — g (x?+,t)dt, j = 0, 1, . . . , n, (4.5) 

then / is a single biquadratic polynomial in every cell of 3 ' . 

Otherwise, suppose in some cell Cj, / is piecewise with horizontal discontinuous-lines y = y^, . . . ,y = 
yi n and vertical discontinuous-lines x = , . . . , x = Xi m of the partial derivatives of order two. Let 

= UiiVh ) • • • j Vim }' ^ = \Ji{ x h >■■■■> x i m , }■ Tnen according to the assumption, c € x \J % is non- 
empty. Without loss of generality, we assume ^ is non-empty. Let y = min y. Suppose Cj is the leftmost 

V&Vy 

cell which takes y = y as its inner discontinuous-line of the partial derivatives of order two with respect to 
x. Denote the x coordinate of the left boundary edge of Cj is Xq, and (xi , y) is a point on the discontinuous- 
line, where xi > xq. Then it follows that 

rx q f Xo d _ 

/ ^-9{s,y-)ds = —g(s,y+)ds. 
J-oo oy J^oo dy 

Since (x\, y) is a discontinuous point, it follows that 

f Xl d f Xl d 

/ Tr9(s,y-)ds ^ / —g(s,y+)ds. 
J-oo oy J-oo dy 

Hence 

d f Xl d 

—g(s,y-)ds^J —g(s,y+)ds. 

This contradicts with the fact that -^-g(x, y) is continuous inside the cell c*. Therefore we have proved that 
/ is a single biquadratic polynomial in every cell of 2T . 

In order to prove / 6 S(2, 2, 1, 1, one needs further to verify that / satisfies HBC. In fact, with 
the same approach, one can prove that / is a single biquadratic polynomial outside 3T. On the other 
hand, f(x, y) is zero when x ^ xq , y ^ yo, where (xq, yo) is the coordinate of the left-bottom corner of 
the T-mesh 3" . Hence / is everywhere zero outside ST , which ensures that / satisfies HBC. Therefore 
/GS(2,2,1,1,^). ■ 

According to the former lemma and the fact dim S(l, 1, 0, 0, S?) — V + , in order to ensure / £ 
S(2, 2, 1, 1, S?\ there are E+A constraints with V + under-determining coefficients, where E is the number 
of interior l-edges. But these constraints are not linear independent, as stated in the following lemma. 

Lemma 4.2 Given a regular T-mesh 3, whose occupied rectangle is (xi, x r ) x [y^, y t ). Assume the dif- 
ferent y coordinates from all the horizontal l-edges are yo < yi < ■ ■ ■ < y n - For g G S(l, 1, 0, 0, 3), it 
follows that 

rx r g rxr g 

J ]fy~ 9 ( s > yt ^ ds = J g^9{s,Vi+)ds, 1 = 0, (4.6) 
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is equivalent with 

g(s,yi)ds = Q, i = l,...,n-l. (4.7) 



77;e similar conclusion can be made for vertical l-edges. 

Proof: We first prove the necessity. Getting started from the bottom boundary 1-edge Iq we have 

■z-9(8, yo-)ds = 0, / g(s, y )ds = 0. 
x, dy J Xl 

Then according to Equation ( 14.61 ) as i = 0, one gets 

d 



On the other hand, according to the piecewise bilinear definition of g, one has 

-9(s, y +)ds = — ff(s, 2/0) H 5(s, yi) 

/ Jx, oy\yi-y y x - y J 

1 



x, dy- 



Vi - Vo \ J Xl 
1 



ds 

(s,yi)ds— I g(s,y )ds) (4.8) 



yi - yo J x 



g(s,yi)ds. 



Therefore 



g{s,yi)ds = 0. 



Recursively, one can prove that Equation ( 14.7b holds for every horizontal 1-edge. 

In order to prove the sufficiency, we take the similar deduction in Equation ( 14.8b . One obtains 

^-g(s,Vi-)ds = ( / g(s,yi)ds- [ g{s, yi-\)ds J , 



x, "y yi - yi-i 

d -g{s,y l +)ds = ( / g(s,y i+1 )ds- I g(s,yi)ds 



Xl 



dy yi+i-yi\Jx l 

Therefore Equation ( 14.71 ) ensures Equation ( 14.61 ) holds. 

Remark: In fact, from the proof of the lemma, we can also conclude that, for horizontal l-edges, 



is equivalent with 



nXr Q rx r Q 

J Yy g(s, Vi -)ds = J —g(s,yi+)ds, i = 0, . 

rx r Q 

J Q^3(s,yi-)ds = 0, i=l,...,n-l. 



Lemma [4721 states that there are at least two redundant constraints among those corresponding to hori- 
zontal l-edges and vertical l-edges, respectively. The following lemma tells us that we can take the constraint 
along four boundary l-edges as redundant ones. 

Lemma 4.3 Given a regular T-mesh ST, the rectangle occupied by it is (xi, x r ) x (yb, yi), where the differ- 
ent y coordinates from horizontal l-edges are yo < yi < ■ ■ ■ < y n - Let g 6 S(l, 1, 0, 0, 2F\ Then Equation 
( 14.61 ) is equivalent with 

d f Xr d 

-K-g(s,yi-)ds = —g(s,yi+)ds, 1 = 1, . . . , n - 1. (4.9) 

Xl dy J Xl dy 

The similar conclusion can be made for vertical l-edges. 
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Proof: The necessity is obvious. In order to prove the sufficiency, we show that Equation ( 14. 9t implies that 
Equation ( 14.71 ) holds. In fact, suppose 

pX r 

g(s,Vi)ds,i = 0,1, ...,n. 

Then Jo = I n = 0. Our object is to prove that, for any i = 1, 2, . . . , n — 1, li = 0. Taking the similar 
deduction in Equation ( 14.81 ), one gets 

t \J Vi+l-Vi f Xr i s , Vi-Vi-1 f Xr / \j 
g(s,y r )ds = / y»-i) H / g(s,y i+1 )ds, 

Vi+i - Vi-i Jx, Vi+i ~ Vi-i Jx, 



i.e., 



T j/i+l ~ j/t r . Vi-Vi-l j 
±1 — ii— 1 

Vi+i - Vi-i Vi+i - Vi-i 



holds for any i = 1, 2, . . . , n — 1. Hence the point set {(t/j, /i)}T=o i s collinear, which means that there 
exists a linear function f(y) such that /(j/i) = /i, i = 0, . . . , n. Since yo — y n = 0, it follows that 
/(y) = 0, which means that It = 0, i = 1, 2, . . . , n. Then the sufficiency is proved. I 



4.2 A lower bound of dimensions 

For a given regular T-mesh £F, which has E interior 1-edges, we assume its associated tensor-product mesh 
£? c has E' interior 1-edges, where E 1 ^ E. Suppose the horizontal 1-edges l\, . . . , in 3T lie on the 
horizontal 1-edge I in 3~ c , where the vertical coordinate of / is y. Then one gets 

/ Ify 9 ^' y~^ ds = / §y~ 9< " S ' ^ + ) ds ' i = lj • ' • ' k 

is equivalent with 

f d f d f d f d 

J dy 9 ^' V ~^ ds = J ~Q^9{ s iV+) ds , i = l,...,k-l, J —g(s,y-)ds = J —g{s,y+)ds. 

According to Lemma [4T2l the constraints, along all the interior horizontal 1-edges and top/bottom boundary 
1-edges in <5F C , 

is equivalent with the integrations of g along all the interior horizontal 1-edges in 3 rc are zero: 

g(s,y)ds = 0. 

The similar conclusion can be made for vertical 1-edges. Therefore the number of the sufficient and neces- 
sary constraints that ensure 1(g) G S(2, 2, 1, 1, 5") for g e S(l, 1, 0, 0, 3f) is just E. 

Moreover, in the tensor-product mesh S? c , if the horizontal knots are xq < x\ < ■ ■ ■ < x m , and the 
vertical knots are yo < Vi < • • ■ < Vm tnen we nave 

/+oo r + oo m-\ n-l 

j g(s,t)dsdt = ^2 _ x i-x) E i = y^(?/j+i - Vj-i)Fj, (4.10) 

where 

Ei = I g(xi,y)dy, F 3 = I g(x,y 3 )dx. 
Jyo J x 

Hence, among the former E constraints, the element number in any maximal linearly independent subset is 
at most E — 1 when these constraints are not empty (i.e., V + > 0). 

Based on these analysis, a lower bound of dimensions of spline spaces S(2, 2, 1, 1, can be proposed 
as follows: 
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Theorem 4.4 Given a regular T-mesh with V + > crossing vertices and E interior l-edges, it follows 
that 

dim S(2, 2, 1, 1, S) > V + - E + 1. 

Proof: Since V + > 0, the constraints are not empty, and E > 1. According the former analysis and 
Theorem l3.41 the dimension of S ( 1 . 1.0,0. 3F) is V + . For any g £ S(l, 1, 0, 0, ST), the constraints ensuring 
1(g) S S(2, 2, 1, 1, have maximal linearly dependent subsets with element number at most E — l. Here 
both V and 1 are linearly injective. Therefore dim S(2, 2,1,1, ,T) ^ V+ - E + 1. ■ 

The lower bound in the theorem can be reached for spline spaces over some T-meshes. For example, 
consider the T-mesh 3T as shown in Figure|2] The B-net method in [3 1 can show that dim S(2, 2, 1, 1, 5T) = 
1. On the other hand, in 3 ', V + — 6, E = 6. Then V + — E + 1 = 1, which means that the lower bound is 
reached in this T-mesh. 

Furthermore, it is easy to verify that, for tensor-product meshes, the former lower bound is exactly the 
same as dimensions of biquadratic spline spaces over the meshes, if V + — E + 1 ^ 0. In the next section, 
we will prove that, for hierarchical T-meshes, the lower bound can be reached in some cases as well. 

5 Dimensions of Spline Spaces S(2, 2, 1, 1, over Hierarchical T- 
meshes 

In this section, a careful analysis on the constraints in Section [4] will help us to build a dimension formula 
of the spline space S(2, 2, 1, 1, ST) over a hierarchical T-mesh. Here the key procedure consists of the 
following components: 

1 . A general hierarchical T-mesh is divided into some so-called crossing-vertex-connected branches 
(Definition l5.U . Then the spline space over the hierarchical T-mesh can be divided into the direct sum 
of some subspaces, each of which is defined over a crossing-vertex-connected hierarchical branch, 
which is a T-mesh as well. See Subsection l5.6l 

2. In a crossing-vertex-connected hierarchical T-mesh, the constraint set is proved to with defective rank 
exact one by the following processing: 

(a) The constraints are converted into a new form to reflect the level structure of the hierarchical 
T-mesh. See Subsection l5.2l 

(b) A new set of basis functions of S(l, 1, 0, 0, Sf) is defined according the structure of the T-mesh, 
such that the occurrence of the basis function coefficients in the constraints is regularized. See 
Subsection l5.3l and Proposition |5.9l 

(c) All the l-edges and the corresponding constraints are ordered according to the structure of the 
T-mesh. For each of the l-edges or the constraints, a characteristic vertex is introduced. Write 
all the constraints in a vector form (Co, C\, . . . , Ct) t in the increasing order. The constraint 
Co can be removed since we have known all the constraints are with defective reank at least 
one. Hence we need to shown that (C\, . . . , Ct) t is with full rank. See Subsection l5.1.3l 

(d) Assume the characteristic vertex of Cj is Vi, i = 1, . . . , T. Arrange the coefficients into a vector 
(Pi, . . . , (3t, @t+i, ■ ■ ■ , Pm) T such that is the coefficient of the basis function, corresponding 
toVi, i = 1,2,..., T, in S(l, 1,0,0,^). 

(e) Then there exists a matrix M = (m^ )txm such that (C%, . . . , Ct) T — M(/?i, . . . , (3m) T ■ It 
can be shown that my = 0, i > j, and ma ^ for any i and j. 

Hence the set of all the constraints is with defective rank one. 
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a. Level b. Level 1 c. Level 2 

Figure 6: A hierarchical T-mesh 



5.1 Hierarchical T-meshes 

A hierarchical T-mesh J4) is a special type of T-mesh which has a natural level structure. It is defined in a 
recursive fashion. One generally starts from a tensor-product mesh (level 0). From level k to level k + 1, 
one subdivide a cell at level k into four subcells which are cells at level k + 1. For simplicity, we subdivide 
each cell by connecting the middle points of the opposite edges with two straight lines. Figure|6]illustrates 
the process of generating a hierarchical T-mesh. For a hierarchical T-mesh , in order to emphasis its level 
structure in some cases, we denote the T-mesh of level k to be 3? k . 

For a given hierarchical T-mesh ST, we can extend it to obtain an extended T-mesh 3* e associated with 
the spline space S(2, 2, 1, 1, 3~\ In the following text, hierarchical T-meshes refer to both the classical 
hierarchical T-meshes and their extension. 

Over a hierarchical T-mesh 3", the dimension of S(2, 2, 1, 1, 3?) may be greater than the lower bound 
in Theorem l4.4l For example, consider the hierarchical T-mesh as shown in Figure [7] where the mesh of 
level is a tensor-product mesh with size 3x3, and in the mesh of level 1, there exists only one cell that is 
subdivided. In this mesh, V + = 5, E = 6. The lower bound is V + — E + 1 = 0. But the dimension of the 
biquadratic spline space over the mesh is at least one obviously. 



Figure 7: A hierarchical T-mesh where the dimension is greater than the lower bound 



5.1.1 Crossing-vertex connected 

In order to ensure the dimensions of spline spaces over hierarchical T-meshes reach the former lower bound, 
we need to focus on a special type of hierarchical T-meshes. 

Definition 5.1 For a regular T-mesh, if, between any two different crossing vertices, there exists a continu- 
ous poly-line, which consists of edges in the mesh, such that every joint between two neighboring horizontal 



17 



and vertical edges on the poly-line is a crossing vertex in the mesh, then the T-mesh is called crossing-vertex 
connected. Such poly-lines are called paths between the crossing-vertices. 
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Figure 8: A hierarchical T-mesh with a path between two connected crossing vertices 



For example, two crossing vertices, labeled with light-gray dots, are selected in the T-mesh as shown in 
Figure[8] A path between these two vertices is illustrated in light-gray as well, where the joints are shown 
with light-gray triangles. 

Definition 5.2 From level k to level k + 1 as forming a hierarchical T-mesh, if there exists a cell of level k 
to be subdivided, but all its horizontal and vertical neighboring cells of level k remain unchanged, then the 
cell is called an isolated subdivided cell. 



For a hierarchical T-mesh, it is crossing-vertex connected if and only if, in any level of forming the hierar- 
chical T-mesh, there do not exist any isolated subdivided cells. 

In Subsections I5.2J - I531 we will prove that the lower bound of dimension proposed in Theorem 14.41 is 
exactly the dimension of S(2, 2, 1, 1, 3?) over a crossing-vertex connected hierarchical T-mesh. Then in 
Subsection 15. 61 a dimension theorem can be proposed to calculate the dimension of the spline space over 
a general hierarchical T-mesh by dividing it into the union of some crossing-vertex connected hierarchical 
branches. 



5.1.2 Level numbers of edges, 1-edges, and crossing-vertices 



edge of level pi — 
edge of level po 
edge of level pi — 















Vl 


V2 













edge of edge of 
level po level pi 



Figure 9: Level numbers of edges, 1-edges, and crossing vertices 



We assign every edge or l-edge in a hierarchical T-mesh with a level number as same as the level number 
of the T-mesh where the edge or l-edge just appears. An l-edge of level k consists of edges of level k. The 
extension of an edge in the extended T-mesh is assigned with the same level as its source. The new adding 
1-edges in the extended T-mesh are assigned with level 0. 
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A crossing vertex is assigned with two level numbers, denoted as (kh, k v ), corresponding to the level 
number kh of the horizontal 1-edge and the level number k v of the vertical 1-edge where the vertex lies, 
respectively, kh and k v are called horizontal level and vertical level, respectively, of the crossing vertex. 

For example, in the hierarchical T-mesh & as shown in Figure [9J suppose the middle horizontal 1-edge 
is with level po, and its upper and lower neighboring 1-edges are with levels p\ and p2, respectively. Here it 
should be po >pi,po > p2- Because v\ is the intersection of two 1-edges with level po, its level is (po,po). 
V2 is the intersection of a horizontal 1-edge of level po and a vertical 1-edge of level p\. Hence its level is 

(po,Pi). 



5.1.3 An ordering on interior 1-edges 

In order to sort the constraints reasonably, we introduce a partial ordering on interior 1-edges in a hierarchical 
T-mesh. This order will be used in the proof of Theorem |5.10| to facilitate the rank analysis of the constraints. 
Before that, we propose the following definition. 

Definition 5.3 In a hierarchical T-mesh 2? , two interior l-edges are continuous if they intersect in a cross- 
ing vertex of ST. An l-edge set S is connected if for any two l-edges £q and £\ in S, there exists a continuous 
series of l-edges lo,h, ■ ■ ■ ,lk in S between £q and l\, i.e., Iq — £q, Ik — l\, and U and li + \ are continuous 
for any i = 0, 1, . . . , k — 1. 
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Figure 10: A hierarchical T-mesh with its decomposition according to the edge levels 



Consider a crossing-vertex connected hierarchical T-mesh ST , Fix a level number fc 0. Then all the 
1-edges with level k is possibly not connected. See Figure [10] for an example, where all the edges of level 
2 are not connected. We assume that they form some maximal connected branches. Denote these branches 
to be 8k.i, i = 1, . . . , Ifc. Hence in the example of Figure [Tol we have To = 1, T\ = 1, and T2 = 2. 
Because 8? is crossing-vertex connected, when k > 1, there exists at least one crossing vertex with level 
number (fc, j) or (j, k) on some 1-edge in the branch 9k.i, where j < k. See FigurefTOlfor examples, where 
the specified crossing vertex is shown in light-gray. This crossing vertex is the intersection between two 
1-edges of level k and j, where the 1-edge of level j is called an entering 1-edge of the branch, which is 
connected with every 1-edge in (9^. The entering 1-edge of all the 1-edges of level zero is defined to be any 
1-edge of level zero. 

Fixed an entering 1-edge £q of 9k,i- For an 1-edge £ in 9k.i, there exists many continuous series of 1-edges 
in 9k.i connecting £ and £q. The number of the 1-edges in a series is called the length of the series, and the 
minimal length of all the series connecting £ and £q is called the distance between £ and £0, denoted 
as dist(^,^o)- Suppose eo = Iq, e±, . .. , e s = £ is an 1-edge series with the minimal length among the 
continuous series connecting £ and £q. Then the intersection point between e s = £ and e s _i is a crossing 
vertex on £, which is called a characteristic vertex of the 1-edge £. Its level is (fc, j) or (j, k), where j ^ k. 
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After having selected the entering 1-edges for all the connected branches 9k,i, we introduce a partial 
ordering <i on all the interior 1-edges in 2? . For any two interior 1-edges l\ and 1% with levels k\ and fo, 
respectively, where £j is in the branch i = 1,2, we define l\ <x £2 if 

1. ki < k 2 , or 

2. fci = fc 2 and i\ < or 

3. ki = k2 and i\ = %2, dist(£i,£o) < dist(£2,^Q), where ^ is the entering 1-edge of the connected 
branch in which lx and £ 2 lies, since these two 1-edges are in the same connected branch. 

This order is not total, because, in Case 3, it is possible that dist(£i,^o) = dist(^ 2 ,^o) f° r two different 
interior 1-edges. 



5.2 Conversion of constraints 

In Section|4]we have proposed some different versions of the necessary and sufficient conditions that ensure 
1(g) € S(2, 2, 1, 1, &) for any g E S(l, 1, 0, 0, ST). In order to facilitate the latter analysis, we take the 
following notation and the constraints. 

For any regular T-mesh ST, denote its occupied rectangle to be (xi,x r ) x (y^, y t ). is an associated 
tensor-product mesh with & . Assume the y coordinates of all the horizontal 1-edges in 3F C to be y^ = y < 
yi < ■ ■ ■ < y m < y m +i — Vt> an d the x coordinates of all the vertical 1-edges in 3~ c to be X[ = x n < X\ < 
■ ■ ■ < x n < x n+ i — x r . For i — 0, 1, . . . , m + 1, assume the horizontal 1-edges . . . , l^ a . are with y 
coordinate?/ = j/,. HereZ|\, . . . ,Zf Q . are sorted from left to right. On the other hand, for j =0,1,... ,n+l, 
assume the vertical 1-edges Vj X , . . . , l^p. are with x coordinate x = Xj. Here . . . , l^ j3 _ are sorted from 
bottom to top. It follows that o.q = a m+ i = [3o = /3 n +i = 1, cti H h a m + (3\ + ■ ■ ■ + (3 n = E. 

With these notation, according to Lemma POl we can allocate the E constraints to the interior 1-edges 
in 3~ c in the following fashion: To any interior horizontal 1-edge y = y,-, the corresponding constraints are 

d f d 

-^-g(s,Vi-)ds = / —g(s,yi+)ds, k = 1, ...,0;,. (5.1) 

To any interior vertical 1-edge x — xj, the corresponding constraints are 

d f d 

Furthermore, by applying the similar deduction with Equation (14.8b , Equations ( I5.lt and ( 15.2b are equivalent 
with 

r x iki r*iki r x iki 

(Vi+x-Vi-i) g(s,Vi)ds = (yi+i-Vi) g(s,yi-x)ds+ (yi-yi-i) g(s,y i+ i)ds, 



k = 1, . . . ,a if (5.3) 

rt/jki rVjki rVjki 

(xj+i—Xj-i) g(xj,t)dt = (xj+i - Xj) g(xj-i,t)dt+(xj —Xj-i) g(x j+1 ,t)dt, 

Jy,ko "JVjho Jvjko 

k = l,...,p j , (5.4) 

respectively, where x^ k0 and x\ kl are the x coordinates of the two end-points of l^ k , and yj fc0 and y'j kl are 
the y coordinates of the two end-points of lj k . 

Now we focus on hierarchical T-meshes, on which the corresponding constraints can be converted into 
a form to reflect the level structure of the T-mesh. At first, we introduce two definitions. 
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In a hierarchical T-mesh select any horizontal 1-edge I with level k > 0. Hence / appears in & since 
3^ k . On I, there exists one or more crossing vertices with vertical level k. These crossing vertices are the 
center of some inserted crossing from 2 rk ~ 1 to 3 rk . The 1-edge I consists of the horizontal edges of these 
inserted crossing. It follows that the vertical edges of the inserted crossing intersect with two 1-edges /' and 
l b in T k ~ x . Here we assume l b with level k b lies under I, and i* with level fc* lies above I. 

Definition 5.4 For a horizontal l-edge I with level k > 0, the l-edges l b and /' are defined as the above 
description. Then l b and l l are called the support l-edges of I. For an horizontal l-edge I with level 0, its 
support l-edges are defined to be the two nearest horizontal l-edges with level 0, which lie above and under 
I, respectively. For vertical l-edges, the support l-edges are defined in a similar way. 

It is obvious that, for any horizontal/vertical 1-edge I, two vertical/horizontal l-edges, which are through 
the two end-points of I, intersect with the two support l-edges of I, where there is not other crossing vertex 
between the intersection points and the end-points along the two vertical/horizontal l-edges. 

Definition 5.5 For the E constraints ensuring 2(g) £ S(2, 2, 1, 1, the number r of its maximal linearly 
independent subset is called the rank of the constraints. Here E — r is called the defective rank of the 
constraints. 

In fact, if we specify a set of basis functions of S(l, 1, 0, 0, and undetermine the coefficients of g 
with these basis functions, then these constraints can be represented into the linear combinations of these 
coefficients. Hence the rank of the constraints is the rank of the corresponding coefficient matrix. 

In Equations ( 15.31 ) and ( 15.4b . a constraint along an 1-edge is represented into the linear combination 
of three integrations, which are along the current 1-edge and its two neighboring horizontal/vertical lines, 
respectively, with the same integration limits. Here the horizontal/vertical lines share the same ylx coordi- 
nates as the two-sided nearest horizontal/vertical l-edges to the current 1-edge. In the following lemma, we 
will convert these constraint into a new form, such that every constraint along an 1-edge is a linear combina- 
tion of three integration, which are along the current 1-edge and its support l-edges, respectively. With this 
form, every constraint will involve undetermined coefficients in a way easy for rank determining. 

Lemma 5.6 Given a hierarchical T-mesh 27 , select an interior horizontal l-edge if. Suppose the y coordi- 
nate of if is y\, the x coordinates of its two end-points are x\ and x[ l 2 , and the y coordinates of its support 
l-edges are yf b and yf*. Then Equation (15 ,3b holds for all the interior horizontal l-edges if and only if, for 
all the interior horizontal l-edges, the following equation holds: 



(vT - Vf) \ 9(8, yf)ds = (yf - yf) / g(s, y[ lb )d s + [y\ - yf) / g{s, yf)ds. (5.5) 



For interior vertical l-edges, the conclusion is similar. 

Proof: Necessity. Suppose the interior horizontal 1-edge l h is with level k, and its two support l-edges are 
l b and /*. Between l h and l b , there are some other l-edges. Suppose the different y coordinates these l-edges 
are yo < < • ■ ■ < y~n, where yo and y fi correspond to l b and l h , respectively. Suppose again the different 
y coordinates of the l-edges between l h and P are y n < yn+i < ■ ■ ■ < yn+m, where y~n+m corresponds 
to I*. According to the definition of the support l-edges, these horizontal l-edges excluding £ , l b and l f are 
with level greater than k. Hence there does not exist one of them such that its vertical projection onto l h 
takes one of the end-points of l h as its interior point. Therefore, according to the constraints in Equation 
(15.31 l corresponding to these l-edges, we can conclude that (j/i-i, h-i), (y~i, h), (§i+x, li+i) are collinear, 
where i = 1, . . . , m + n — 1, and 







21 



Then (j/q, Io), (yn, In), (lJh+m, In+m) are collinear, which means the corresponding constraint in Equation 
Ol holds. 

Sufficiency. We prove it in an inductive fashion from the highest level to the lowest level. We will show 
that, for any given level ko, if the constraints as defined in Equation ( 15.51 ) hold corresponding to the 1-edges 
of level k ^ fco, then the constraints in Equation ( 15.3b corresponding to all the 1-edges of level fco hold as 
well. 

Suppose the maximal level number in the hierarchical T-mesh £T is M. For an interior horizontal 1-edge 
with level M, since there do not exist other 1-edges between the current 1-edge and its support 1-edges, the 
corresponding constraint defined in Equation d5.3t is the same as one defined in Equation ( 15.51 ). 

Now we assume that, for all the interior horizontal 1-edges with level greater than fc, if the corresponding 
constraints in Equation ( 15.51 ) hold, then the corresponding constraints in Equation ( 15.3b hold as well. Select 
an arbitrary interior horizontal 1-edge l h with level k, whose support 1-edges are l b and I*. The two-sided 
nearest 1-edges of l h are with y coordinates ho and hi. Since the level numbers of the other horizontal 
1-edges between l h and l b are greater than k, according to the inductive assumption, the integration values 
along these horizontal lines with the integration interval , x^ 2 ] are collinear with respect to their y coor- 
dinates. Especially, the corresponding integration values along i , l b and y — ho are collinear with respect 
to their y coordinates. Similarly, the integration values along l h , l l and y — hi are collinear with respect to 
their y coordinates as well. Hence the integration values along l h , y = ho and y — hi are collinear with 
respect to their y coordinates. This finishes the proof of the sufficiency. B 

In the next subsection we will specify a proper set of basis functions of S(l, 1, 0, 0, such that the 
coefficients of g under these basis functions appear in these constraints in a regular form. 

5.3 Hierarchical basis functions of S ( 1 , 1 , , , ^) 

In Subsection l3.4l an approach has been proposed to specify a set of basis functions of S(l, 1, 0, 0, 3") with 
many good properties. Now we specify a new set of basis functions for S(l, 1, 0, 0, ST), which does not 
poss the unity -partition property. But under this set, the former constraints in Lemma IBT61 will appear in a 
regular form, which facilitates the rank determination of the constraints. 

The new set of basis functions is defined level by level when forming the hierarchical T-mesh. Every 
basis function is associated with a crossing vertex. At the level 0, we consider the level T-mesh 
and introduce the functions to be standard bilinear tensor-product B-splines. Hence every function can be 
associated with a crossing vertex in in a way that the function reaches one at the crossing vertex, and 
zero at all the other crossing vertices of 3?° . Use the set B° to denote all these functions. Suppose the 
current level number is k ^ 1. Consider a new-coming crossing vertex in this level. 

1 1 1 1 1 1 

i i 

( , i I i 

I 1 1 

i i 

I I I I I I 

a) A part in 3T k ~ x b) After inserting crosses c) Support example 1 d) Support example 2 

Figure 1 1 : The support of the hierarchical basis functions 

1. If its level is (k, k), then the crossing vertex must be the center of an inserted cross into a cell c of 
T . Hence a function can be defined associated with the crossing vertex such that the function 
reaches one at the vertex and its support is c. See Figure[TT]c) for an example, where the filled region 
with light-gray is the support of the specified function associated with the new crossing vertex labeled 
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with •. Here the discontinuous of the derivatives d/dx and d/dy appears on the edges of the inserted 
cross. 

2. Otherwise, the crossing vertex is the middle point of an edge e in 3 rk ^ 1 , where e is the common edge 
of two neighboring cells c\ and C2, each of which is subdivided by inserting a cross from level k — 1 
to k. Then a function can be defined such that the function reaches one at the current crossing vertex, 
its support is exactly c\ U C2, and its derivative discontinuous in the support lies only on the lines 
through the current crossing vertex. See FigureQT|d) f° r an example, where the derivatives d/dx and 
d/dy are continuous along the dashed edges. 

All the functions introduced at level k are denoted to be B k . 

Lemma 5.7 Suppose the maximal level number of the hierarchical T-mesh S? is M. Define B = (Jfc=o 
Then B is a set of basis functions o/S(l, 1, 0, 0, S?\ 

Proof: According to the definition of the functions in B k , it follows that the number of the functions 
in B is exactly the number of crossing vertices in 27 " . Moreover, these functions are in the spline space 
S(l, 1, 0, 0, Therefore, in order to prove that B forms a set of basis functions of S(l, 1, 0, 0, 3^), one 
just needs to show the functions in B are linearly independent. 

Assume the region occupied by 3*" is £1. Let B k = {b k (x, y), . . . , b\ (x, y)}. Suppose a set of coeffi- 
cients of a k ensure that 

M ri fe 

fix, y ) :=EE a * b *( x > y) = °< (*> y) e n - 

fc=0 i=0 

Now we will prove that a k = for any i and k. 

Consider the values of / at every crossing vertex vf with level (0, 0). Since / = 0, it follows that 
f( v i) = 0- O n tne other hand, the function bf(x : y) in £>° associated with is one at v®, and all the other 
functions vanish at v®. Hence /(«§) = a i^i = i.e., a® = 0. Hence, we have 

M n k 
k=l i=0 

Suppose we have proved that a\ = 0, i = 0, . . . , nj, j = 0, 1, . . . , k— 1. Now we consider the functions 
and their coefficients introduced at level k. Each of such the functions is associated with a crossing vertex. 
We first consider the functions which are associated with crossing vertices whose vertical and horizontal 
level numbers are different. On such a crossing vertex, all the functions are zero except the associated 
function with the current crossing vertex. Hence the corresponding coefficient is zero. Then we consider 
the other functions, which are associated with crossing vertices with the same vertical and horizontal level 
numbers. One can conclude that their coefficients are zero as well in a similar way. Hence all the coefficients 
are zero, and the functions in S3 are linearly independent. Therefore they form a set of basis functions of 
£3(1,1,0,0, BT). ■ 

W will call this set of basis functions to be the hierarchical basis functions of S(l, 1, 0, 0, 3 r ). 
5.4 Concretion of the constraints 

For a given hierarchical T-mesh 2T, suppose the hierarchical basis functions of S(l, 1, 0, 0, 5*) are 

{b>i(x,y)} = {cj&y)}^, 
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where b!- (x, y) is associated with a crossing vertex appearing in & since level k. 
Represent g as 

v+ 

g{x,y) = ^2ajCj(x,y). (5.6) 

3=1 

Since there is a one-to-one mapping between the basis functions and the crossing vertices, the coefficient 
of a basis function is called also the coefficient of the corresponding crossing vertex. 

Now we state the characteristic of the constraints after substituting the representation of g as defined 
in Equation ( 15.6b into Equation (15.5b which ensures 1(g) € S(2, 2, 1, 1, At first we consider the 
constraints corresponding to the interior 1-edges of level 0. 

Proposition 5.8 After substituting Equation d5.61 > into Equation d5.5l ), and then taking a proper transfor- 
mation, all the constraints corresponding to the interior horizontal l-edges of level are in such a form that 
the nonzero terms in a constraint consist of those associated with the crossing vertices on the same l-edge. 

Proof: Each of interior horizontal 1-edges with level must traverse the whole mesh. Its support 1-edges 
are with level as well. Over the T-mesh we apply Lemma |4~2l and obtain that all the constraints 
corresponding to the interior horizontal 1-edges of level are equivalent with that the integration of g along 
each interior horizontal l-edge is zero. For every l-edge of level 0, among all the coefficients, only those 
of the crossing vertices on the current l-edge are nonzero. Hence the integration can be represented into a 
linear combination of the coefficients of the crossing vertices on the current l-edge. H 

Then we consider the constraints corresponding to the interior 1-edges of level greater than 0. 

Proposition 5.9 Suppose the current interior horizontal l-edge is with level k > 0. After substituting 
Equation (15.6b into Equation (15.5b , the nonzero coefficients of the basis functions in the corresponding 
constraint consist of two parts as follows: 

A. the coefficients of the crossing vertices on the current l-edge; 

B. the possible coefficients of the crossing vertices with horizontal level less than k and vertical level 
greater than k. 

Proof: Select an interior horizontal l-edge l h with level k > 0. Suppose its support 1-edges are l b and i* 
with levels k\ and k%, respectively. Assume the x coordinates of the two end-points of l h are Xq and x\, 
and the vertical 1-edges through the two end-points are I? and l v r . Then the two vertical 1-edges If and 1% 
intersect with 4 and l t . Consider the crossing vertices whose associating basis functions are nonzero on l h , 
l b \[x h x h ] 0T ^\[x h x h }- These crossing vertices can be classified into the following cases with respect to their 
level 

1. k'j < k: 

(a) k\ ^ k. Because, for any x G [xq,Xi], the corresponding three points on i , l f , and l b with 
the horizontal coordinates x appear simultaneously in a cell of 3~ k ~ x (including its boundary), 
it follows that for the current basis function b(x, y), Equation (15.5b holds as g = b(x, y). This 
means that the coefficients of b(x, y) does not appear in the constraint after simplification. 

(b) k\ > k. This types of coefficients can appear in the constraint. 

2. kg — k: This type of crossing vertex has no contribution to the constraints, unless it lies on l h . 

3. fcg > k: The associating basis function with this type of crossing vertex vanishes on l h , l b and Z* 
between Xq and x\. Hence the corresponding coefficient does not appear in the constraint. 
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Therefore the nonzero coefficients come in Cases 1(b) and 2, which correspond to Cases B and A in the 
proposition description, respectively. I 

We can make a similar classification on the coefficients' appearance in the constraints corresponding to 
interior vertical 1-edges. 

After these classifications, one knows that the coefficient of a crossing vertex with level (fci, fo) can 
appear in the following three types of places: 

1 . The constraints corresponding to the horizontal 1-edge through the vertex; 

2. The constraints corresponding to the vertical 1-edge through the vertex; 

3. If ki < &2 — 1 (or k\ — 1 > &2), it mav appear in the constraints corresponding to horizontal (or 
vertical) 1-edges with level less than (or k\). If k\ — k-2 or &2 ± 1, this situation does not happen. 

The involved 1-edges in the first two cases are called naturally appearing 1-edges of the coefficient. The 
involved 1-edges in the third case are called unnaturally appearing 1-edges of the coefficient. For simplify, 
a constraint corresponding to a horizontal/vertical 1-edge of level k is simply called the horizontal/vertical 
constraint of level k in the following. For two given interior 1-edges £\ and £2, their corresponding con- 
straints are c\ and C2, respectively. If £\ <i £2, then we define c\ <\ C2 as well. 

5.5 Dimension theorem over crossing-vertex connected hierarchical T-meshes 

With these preparations, now we state and prove the dimension theorem of biquadratic spline spaces over 
crossing-vertex connected hierarchical T-meshes. 

Theorem 5.10 Over a crossing-vertex connected hierarchical T-mesh S? or its extension associated with 
S(2, 2, 1, 1, ST), where V + > 0, it follows that 

dim S(2, 2, 1, 1, ST) = V + - E + 1. 

Proof: Recall some facts first. For any g G S(l, 1, 0, 0, in order to ensure that T{g) £8(2,2,1,1,5"), 
the constraints corresponding to every interior horizontal 1-edges (defined in Equation ( 15.5V ) and vertical 
1-edges should be satisfied. Here we apply the hierarchical basis functions B of S(l, 1, 0, 0, 3F). These 
constraints can be represented into linear combinations of the coefficients of g under the basis functions B. 
The coefficients appear in these constraints as stated in Propositions 15. 8| and |5. 91 

Suppose the 1-edges of level k in !7 form some connected branch 8k.i, i = 1,2, . . . , Tfc, and the entering 
1-edges for all the branches O^.i have been selected, denoted as Specially, all the interior 1-edges of 
level zero just form a connected branch, and its entering 1-edge is selected to any one interior 1-edge of level 
zero. Then we can introduce a partial ordering <i on the interior 1-edges and the corresponding constraints 
as stated in Subsection l5.1.3l 

At the beginning, all the constraints are linearly dependent because Equation H.lOt shows a linear 
combination of these constraints with result zero. In this linear combination, each of the coefficients of 
level zero is nonzero. Hence the rank of the constraints remains unchanged after deleting any one constraint 
of level zero. Without loss of generality, we assume that the deleted constraint is corresponding to the 
entering 1-edge £0,1 of #0,1 ■ Then we focus on the remaining constraints and we will show that they are 
linear independent. 

Sort all the remaining constraints into a non-decreasing series C\ , C%, ■ ■ ■ , Ct according to the ordering 
<i, where T = E — 1, and sort the coefficients {a^} in Equation ( 15.6b into a series fix, fa, . . . , (3y+ such 
that Pi is a characteristic vertex of the 1-edge whose corresponding constraint is Cj, i = 1, . . . ,T. The rest 
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variables (3t+i, ■ ■ ■ , Pv+ are arranged randomly. Then we can write these constraints into the following 
matrix form: 

(C 1 ,C 2 ,...,C T ) T = M(/3 1 ,/3 2 ,...,/? v+ ) T , 

where M = (m^- ). Because the characteristic vertex of an 1-edge I with level k is with level (k, j) or (j, k), 
where j ^ k, it follows that, according to Propositions l5.8l and |5.91 the coefficient /3j does not appear in the 
constraints C%+\ , . . . , Ct, where i — 1, . . . ,T — 1. Hence = 0, i > j. On the other hand, the matrix 
M is with ma ^ 0, i = 1, . . . , T. This means that the matrix M is with full row rank, i.e., rankM = T. 
Therefore the dimension of the spline space S(2, 2, 1, 1, 3) is V + — T = V + — E+ 1. This completes the 
proof of the theorem. H 



5.6 General Dimension Theorem 

In this subsection we consider the dimension formula of S (2, 2, 1, 1, 3), where 3? is a general hierarchical 
T-mesh. 

At first, we discuss how to divide a general hierarchical T-mesh into the union of some crossing-vertex 
connected hierarchical T-meshes. Suppose the given hierarchical T-mesh is 3F , and let Cj, % = 1, . . . , C, be 
all the isolated subdivided cells with level k > in forming ST , Then the subdivision happening in Cj to 
form 2T will form as well a hierarchical T-mesh, denoted as Here % occupies the same region as the 
cell Cj. In the following we will apply a T-mesh operation '\'. Suppose i = 0, 1, . . . , k, are hierarchical 
T-meshes. Then 3q\{3i, . . . , 5^} is a new T-mesh which consists of edges and vertices in ,% but not in 
the interior of 31, i = 1, . . . , k. Let % — 3~ , and 

% = %\{% : % is a submesh of %, j = 1, . . . , C, j ^ i}, i = 0, 1, . . . , C. 

Then it is easy to verify that % is a crossing-vertex connected hierarchical T-mesh. Here % comes from 
5T, and the other "%'s come from the isolated subdivided cells with level k > 0. Then it follows that 3" can 
be seen as the disjoint union of i = 0, 1, . . . , C. 

Lemma 5.11 Suppose the hierarchical T-mesh S? and its disjoint union of i = 0,1, ... ,C, are defined 
as before. Then it follows that 

c 

8(2,2,1,1,^) =08(2,2,1,1,*). (5.7) 
t=o 

Proof: At first, we prove that the intersection of any two different subspaces among S(2, 2, 1, 1, i = 
0,1, ... ,C are {0}. Therefore we can define the direct sum of these subspaces. Suppose we select the 
submeshes ^ and where io ^ i\. The submeshes % and ^ t are defined as before. If any one 
of a i/i a and a ^i x is not a submesh of the other one, then the regions occupied by these them are disjoint, 
which follows that the intersection of the subspaces over ^ and ^ is just {0}. Otherwise, we assume, 
without loss of generality, that % is a submesh of . Then the region occupied by % is inside a cell of 
. Denote the cell to be c. In "f^ , the zero function is a unique function whose function values and two 
partial derivatives of order one are zero along the boundary of c. Hence it follows that the intersection of 
the subspaces over ^ and ^ is {0} as well. 

It is obvious that 

c 

8(2,2,1,1,^) 3 08(2,2,1,1,^), (5.8) 

i=0 

On the other hand, for any / G S(2, 2, 1, 1, 3), one can construct its component in each subspace as 
follows. For any j = 0,1, ... ,C, we can arrange all the meshes {^i}f =1 , each of which takes ^ as a 
submesh, in an ascending chain as follows: 

% = % C %-x C---C% =%. 
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Since Yq comes from = & by deleting the subdivisions in some isolated subdivided cells, we can 
define a new function Pq/ which meets / with order one along all the edges in Yq. Then / — Pq/ vanishes 
out of those isolated subdivided cells of We define recursively that Fjf is a function in S(2, 2, 1, 1, Yf) 

which meets the function / — Y^k=i ^ k f w ^ or der one along all the edges in Yj. Then / — Y^k=i ^ fe / 
vanishes in Yj except in the isolated subdivided cells of %j . 

According to the definition of Pfe, it follows that / — J^Pfc/ vanishes everywhere in & . Hence / e 
S(2, 2, 1, 1, Yi), which means that 

c 

S(2, 2,1,1^)C0 S(2, 2, 1, 1, Yi). (5.9) 

i=0 

Combining Equations d5 ,8b and (15.9b together, one gets Equation (15.71 >. I 

With this lemma, we have the following theorem about the dimension formula of the spline spaces over 
general hierarchical T-meshes. 

Theorem 5.12 Suppose 2? is a hierarchical T-meshes with (5 — 1 isolated subdivided cells. Then 

dimS(2, 2, 1, 1, ST) = V + - E + 5. 

Proof: Suppose ^ and Yj, j = 0, 1, . . . , C are defined as in the proof of Lemma l5.HI Then it follows 
that C = 5 — 1 and, according to Lemma B.llI 

c 

dim S(2, 2, 1, 1,£T)=^2 dim S(2, 2, 1, 1, %). 

i=0 

Assume in Y. there are crossing vertices and Ei interior 1-edges. Then 

c 

dim S(2, 2, 1, 1, ^) = YSy? ~ E i + !)• 

Since any two different meshes among Vi do not share any common crossing vertices and interior 1-edges, 
it follows that 

c 

^2(V+ - Ei + 1) = V+ - E + C + 1 = V + - E + S. 

Hence we have 

dimS(2, 2, 1, 1, £T) = V + - E + S. 



5.7 Some Notes on Construction of Basis Functions 

After obtaining the dimension formulae of biquadratic spline spaces over hierarchical T-meshes, we natu- 
rally consider how to construct their basis functions with some good properties as stated for bilinear basis 
functions in Section [3T4l 

To do so, we fist need to make clear of the topological meaning of V + — E + S. 

Definition 5.13 Given a hierarchical T-mesh ST, one can construct a graph §f by keeping all the cross- 
ing vertices and the line segments with two endpoints being crossing vertices, and removing all the other 
vertices and the edges in ST. is called the crossing-vertex-relationship graph ( CVR graph for short) of 
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&x CVR graph of J x 



^ 2 CVR graph of ST 2 

Figure 12: Two examples of CVR graphs. 

See Figure [12] for two examples. Here 07\ has an isolated subdivided cell at its right-top part. Hence 
5 = 2. Its CVR graph has two disconnected parts. £7% is an extended T-mesh of 3\ with respect to the spline 
space S(2, 2, 1, 1, S7\). 3^ has no isolated subdivided cell. The corresponding CVR graph is connected. 
The following theorem states the relationship between V + — E + 6 in 3? and the cell number in . 

Theorem 5.14 Given a hierarchical T-mesh SF with V + crossing vertices, E interior l-edges and 6—1 
isolated subdivided cells, suppose there are Fc§ cells in its CVR graph §f. Then it follows that 

V + - E + 6 = Fcg. (5.10) 

Proof: Suppose there are edges in Sf. According to the definition of the CVR graph, it follows that 
there are 6 disconnected parts in . Hence with the Euler formula, one gets 

Fcg - Ey + V + = 6. 

Hence in order to prove Equation ( |5.10t , one just needs to show that 2V + = E + Ecg, In fact, consider 
any 1-edge £j with crossing vertices. Then 1^ generates — 1 edges in Sf . After running through all 
the l-edges in 3F , each of the crossing vertices is met twice. Hence it follows that 2V + = E + E%, which 
finishes the proof of the theorem. H 

Theorem l5 . 1 4l hints us that the basis functions of the biquadratic spline space over a hierarchical T-mesh 
could be constructed around the cells in its CVR graph. Some experiments have been done on this idea, 
which will be explored in the future. 

Furthermore, we expect that CVR graphs will play an important role in dimension analysis and basis 
function construction of higher degree spline spaces over (hierarchical) T-meshes. For example, we have 
the following conjecture: 

Conjecture 1 Suppose 3 is a hierarchical T-mesh, and its CVR graph is Sf. As m,n 2, it follows that 

dim S(m, n, m — 1, n — 1, 2?) = dim S(m — 2, n — 2, m — 3, n — 3, 
where the spline space S(m, n, a, /3, &) is defined in a similar way with the spline space over a T-mesh. 
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Theorem 15 . 1 41 states that the conjecture holds as m = n = 2. As for rn = n = 3, we have tried many 
examples, which support this conjecture as well. 



6 Conclusions 

In this paper the dimension of bilinear and biquadratic spline spaces over T-meshes are discussed. The 
basis strategy is by linear space embedding with an operator of mixed partial derivative. We obtained the 
dimension formula of bilinear spline spaces over general T-meshes, and that of biquadratic spline spaces 
over hierarchical T-meshes. Only a lower bound of the dimension is build for biquadratic spline spaces over 
general T-meshes. 

In the future, the basis function construction of biquadratic splines spaces and the proposed conjecture 
will be explored. 
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